AD-A0B1  179  ROCKWELL  INTERNATIONAL  THOUSAND  OAKS  CA  SCIENCE  CENTER  F/6  20/9 
STUDIES  OF  UPPER  SURFACE  BLOWN  AIRFOIL^  IN  INCOMPRESSIBLE  AND  T— ETC (l 
JAN  80  N  0  MALMUTH»  W  D  MURPHY.  V  SHANKAR  N00014-76-C-0350 

nun  Acctcre-n  VMW.ISTB _ NL _ 


w 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  This  PAGE  fWttmn  0«<  Entered; 


O*  ' 

1 

r—4 
1— 4  ft 
GO  ' 
O 
•sll 

Si 


REPORT  DOCUMENTATION  PAGE  /  /  (/ 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

IiIfiiMMHHHIifl 

F3.  RECIPIENT'S  CATALOG  NUMBER 

4.  wad- Subtitle)  —  ‘  *  "V 

STUDIES  OF  UPPER  SURFACE  BLOWN  AIRFOILS  IN 
^INCOMPRESSIBLE  AND  TRANSONIC  FLOWS 

m  ...  .  *  >>  — 

i  /  Gy ) 

5.  TYPE  OF  REPORT  8  PERIOD  CO  'ERED 

Technical  Report 

09/30/78  through  09/30/79 

Y  AUTHOR/-.)  1  /  J  7  ' 

7  N.D./Malmuth*  W.D. /Murphy/  V.^hankar/j.D./Cole 
(uu.m)  ,~a7Ul^  Ey  Cumberbatch  JPurdue  University) 

8  CONTRACT  OR  GRANT, NUMBER/*; 

N0O914-76-C-O350, 

Vs) 

9.  PERFORMING  ORGANIZATION  NAME  ANO  AOORESS 

Rockwell  International  Science  Center'- 
P.0.  Box  1085  " 

Thousand  Oaks,  California  91360 

19.  PROGRAM  ELEMENT.  PROJECT.  TASK 

AREA  i  *ORK  UNIT  NUMBERS 

Task  061-234 

11.  CONTROLLING  OFFICE  NAME  ANO  AOORESS 

Office  of  Naval  Research  v’  ,rJ  i 

tz.  REPORT  DATE 

Jani— 

800  North  Ouincy  Street  "V  if, 

Arlington,  Virginia  22217 

13.  NUMBER  eF  P  AG'es 

17 

IA.  MONITORING  AGENCY  name  &  AOORESS  (it  dil/erent  Irom  Conlrollina  Oitice) 

_ ! 

15.  SECURITY  CLASS,  (of  this  report) 

UNCLASSIFIED 

15«.  DECLASSIFICATION  DOWNGRADING 
SCHEDULE  _ 

? 


I 


Approved  for  Public  Release;  Distribution  Unlimited 


/  <.<•  •  /»  '  /UA-  f 
r 


?  V  b  < 
^  A  ( 


n  J 


ePV. 


ecte 

>900 


17.  DISTRIBUTION  STATEMENT  /of  the  abstract  entered  in  Block  20,  it  different  from  Report)  9 


JL 


'■TIC  ABLE. 
r.NED  A 

18.  SUPPLEMENTARY  NOTES  ** '•'•‘'u Li.  JALu j  3Jj y. 

AIAA  Paper  AIAA-80-0270,  presented  at  the  AIAA  18th  Aerospace  Sciences 
Meeting,  January  14-16,  1980,  Pasadena,  California 


19.  KEY  WQROS  'Continue  on  reverse  side  if  necessary  and  identify  by  block  number ) 

Upper  Surface  Blowing,  Jet  Flaps,  Propulsive  Lift,  Transoni  .  Flow,  Numerical 
Analysis,  Partial  Differential  Equations,  Programming,  Asymptotic  Expansions, 
Boundary  Layers,  Viscous  Flow 


Ur 


ABSTRACT  /'Continue  on  reverse  side  if  neceseery  end  identify  by  block  number) 

Asymptotic  and  computational  methods  have  been  utilized  to  study  the 
incompressible  and  transonic  flow  over  upper  surface  blown  airfoils.  To 
provide  a  framework  for  more  approximate  simulations  which  are  subsequently 
discussed,  a  full  potential  formulation  is  given  and  various  numerical 
treatments  are  proposed.  In  this  and  the  other  models,  the  problem  has  been 
decomposed  into  the  treatment  of  the  fine  structure  of  the  jet  and  the 
analysis  of  the  flow  outside  of  it.  Asymptotic  expansions  of  limit  process 


y 


dd 


AN  73  1473 


EDITION  OE  I  NOV  8»  IS  OBSOLETE 


UNCLASSIFIED 


_L 


SECURITY  CLASSIFICATION  of  this  PACE  -l*Tien  Dale  Entered' 


.5 


_ UNCLASSIFIED _ 

^SECURITY  CU*SSIFiC*T:oh  THIj  Ptot'iwun  Dtit  Enund) 

^type  have  been  used  to  treat  the  jet  in  a  thin  layer  approximation  using 
suitable  strained  variables.  Although  vorticity  must  be  accounted  for  in 
matching  with  the  external  flow,  its  effect  on  the  Spence  boundary  conditions 
derived  under  irrotational  assumptions  is  nil  in  regions  away  from  the 
trailing  edge  and  jet  exit.  A  similar  conclusion  applies  for  compressibil ity. 
The  condition  of  flow  pressure  and  direction  compatibility  replacing  the 
Kutta  condition  for  the  unblown  configuration  indicates  that  a  dividing 
streamline  leaves  tangent  to  the  upper  surface  of  the  airfoil  at  the  trailing 
edge.  Computational  results  for  a  USB  airfoil  indicate  significant  enhance¬ 
ments  in  lift  with  blowing.  Comparisons  with  experiments  indicate  that 
viscous  wall  jet  effects,  wave  interaction  phenomena  with  the  mixing  zones 
near  the  jet  exit  and  trailing  edge  layers  must  be  incorporated  into  the 
model  for  improved  simulation  of  the  flow  physics. 


1  4  s!  n 


v~,-  »  ,  ; 

■  \  {  -  i 


j_  Sc««al<a 

I  51-  ~~ - 7i/C 

;  etc  ta2 

Qf  ^  i  V  I 

/J«tinc.*ioa  U'' 


LPJs;  •  l; f 

;  ><,„  , 

I  "  t  ■- 

i  j  A va ? i  r:  ,  ... 

|Dfnt  I  fipc-?  al 


Ilk 


UNCLASSIFIED 


$ecu*l*v  Ci.*SSI  *T:CN  3^  *-ll  Sr.itrtd) 


'  DISCLAIMER  NOTICE 

> 


THIS  DOCUMENT  IS  BEST  QUALITY 
PRACTICABLE.  THE  COPY  FURNISHED 
TO  DDC  CONTAINED  A  SIGNIFICANT 
NUMBER  OF  PAGES  WHICH  DO  NOT 
REPRODUCE  LEGIBLY. 


AIAA-  80-0270 


Studies  of  Upper  Surface  Blown  Airfoils  in 
Incompressible  and  Transonic  Flows 


N.D.  Malmuth,  W.D.  Murphy,  and  V.  Shankar 
Rockwell  International  Science  Center, 
Thousand  Oaks,  CA 

J.D.  Cole 

University  of  California,  Los  Angeles,  CA 
and 

E.  Cumberbatch 

Purdue  University,  West  Lafayette,  IN 


AIAA  18th 

AEROSPACE  SCIENCES 
MEETING 

January  14-16,  1980/Pasadena,  California 


For  permission  to  copy  or  republish,  contact  the  American  Institute  of  Aeronautics  and  Astronautics, 
1290  Avenue  of  the  Americas,  New  York,  N.Y.  10019. 


STUDIES  OF  UP»ER  SURFACE  3L0WN  AIRFOILS  III  I'ICCMPRESS ISLE 
AND  TRANSONIC  :LCWSa 


N.D.  Malmuth*,  W.3.  Murony*",  ana  7.  Shankar'* 
Rockwell  International  Science  Center 
'housand  Oaks,  California  91361 

J.3.  Cole* 

University  of  California  at  Los  Angeles 
Los  Angeles,  California  90024 

ana 

£.  CsJnberbatch 
Purdue  University 
West  Lafayette,  Inaiana 


Abstract 

Asymototic  and  ccmoutational  methods  have  been 
utilized  to  study  the  incomoressible  ana  transonic 
flow  over  ^ooer  surface  blown  airfoils.  To  provide 
a  "rafnework  for  more  aooroximate  simulations  which 
are  subseauently  discussed,  a  full  potential  formu¬ 
lation  is  given  and  various  numerical  treatments 
are  orooosed.  In  this  and  the  other  models,  the 
problem  has  been  decomposed  into  the  treatment  of 
the  fine  structure  of  the  jet  and  the  analysis  of 
the  flow  outside  of  it.  Asymototic  expansions  of 
limit  process  type  have  been  jsed  to  treat  the  jet 
in  a  thin  layer  aooroximation  using  suitable 
strained  variables.  Although  /orticity  must  be 
accounted  *‘or  in  matching  with  the  external  flow, 
its  effect  on  the  Soence  boundary  conditions 
derived  under  irrotational  assumptions  is  nil  in 
regions  away  *‘rom  the  trailing  edge  ana  jet  exit. 

A  similar  conclusion  applies  *or  comoressibil ity . 
'he  condition  of  flow  oressure  ana  direction  com¬ 
patibility  replacing  the  <utta  condition  for  the 
unolown  configuration  indicates  that  a  dividing 
streamline  leaves  tangent  to  the  uooer  surface  of 
the  airfoil  at  the  trailing  edge.  Computational 
results  for  a  USB  ai-foil  indicate  significant 
enhancements  in  lift  with  blowing.  Comparisons 
with  experiments  indicate  that  viscous  wall  jet 
effects,  wave  interaction  phenomena  with  the  mixing 
zones  near  the  jet  exit  and  trailing  edge  layers 
must  be  incorporated  into  the  model  for  imoroved 
simulation  of  the  flow  ohysics. 

1.0  Introduction 

Uooer  surface  blowing  (USB)  nas  been  proposed 
as  a  means  of  increasing  usaole  lift  and  thereby 
enhancing  7/STOL  capability  at  low  speeds  in  land¬ 
ing  configurations.  At  transonic  Mach  numbers,  it 
has  the  further  abdication  of  achieving  low  turn 
radii  in  dogfight  scenarios.  The  attendant  high 
accelerations  are  accomplished  threugn  elimination 
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of  separation  by  suporession  of  adverse  pressure 
gradients  in  the  viscous  boundary  layer,  and  also 
movement  of  shocks  downstream  of  the  trailing  eage, 
thereby  discouraging  shock  induced  seoaration  ana 
buffet  at  high  maneuver  incidences.  Further  appli¬ 
cations  of  laminar  flow  control  througn  stabiliza¬ 
tion  using  tangential  blowing  to  achieve  favorable 
pressure  gradients  is  of  strong  interest  currently. 

In  the  application  of  this  conceot,  the 
engine  bleedoff,  thrust,  and  structural  penalties 
required  to  achieve  the  foregoing  aerodynamic 
advantages  is  of  importance  to  the  designer.  To 
obtain  this  relationship,  a  knowledge  of  the 
associated  flow  fields  is  required.  Although 
attention  has  been  given  to  the  jet  flao  in  theo¬ 
retical  Investigations,  -eiatively  little  analysis 
has  been  performed  on  upper  surface  blown  configu¬ 
rations.,  "or  incomoressible  speeds,  the  work  of 
Soence'’ 1  represents  the  classical  thin  airfoil 
t-eatment  of  the  jet  flao  problem.  At  transonic 
Mac.n  numbers,  a  computational  jet  flao  solution 
based  on  small  disturbance  theory  was  develooed  for 
air*oils,  and  generalized  for^ three-dimens ionai 
wings  by  Malmuth  and  '•urshy .  (2-3)  [n  these  analy¬ 
ses,  the  classical  karman  SudeHey  model  was 
aoolied  with  a  generalized  version  of  the  Murman- 
Cole  successive  line  overrelaxation  schemed!  to 
treat  the  *>ee-jet  boundaries.  The  jet  was 
assumed  to  be  thin,  and  it  was  assumed  on  a 
heuristic  basis  that  the  Spence  boundary  conditions 
were  aoolicable  across  it.  These  conditions 
involve  edui 1 ibration  between  tne  normal  pressure 
gradient  ana  the  centrifugal  *‘orce  associated  with 
the  momentum  in  the  jet. 

In  this  paper,  the  aool icaoil ity  of  the  condi¬ 
tions  will  be  analyzed  -‘or  a  compressible  rota¬ 
tional  jet  in  the  context  of  blowing  upstream  of 
the  trailing  edge  on  the  upper  surface,  i.e., 
upper  surface  blowing  in  contrast  to  the  jet  *lao 
conf iguration  in  which  the  jet  emanates  from  the 
trailing  edge,  'urtnermore,  the  thin  airfoil  jet 
flao  problem  formulation  stioulates  an  initial 
angle  of  the  jet  at  the  trailing  edge.  In  actual¬ 
ity,  this  angle  is  a  function  of  tne  jet  exit  con¬ 
ditions  and  the  local  flow  details,  for  jet  'laos 
and  USB.  In  both  cases,  the  conditions  for  the 
trailing  edge  dividing  streamline  represent  a 
generalization  of  the  <utta  condition  for  the 
unblown  case.  The  aspect  of  the  pacer  involving 
fine-structure  of  the  jet  layer  represents  an 
extension  of  the  earlier  work  of  Malmuth  and 
MurphyiSl  an  transonic  wall  jets.  From  these 
analyses,  the  paoer  will  describe  the  numerical 
aoproacn  to  treat  the  US3  problem,  and  various 
results  showing  possibilities  for  lift  augmentation 


will  be  presented.  An  additional  significance  of 
the  work  is  its  relevance  to  the  wake  curvature 
effects  treated  ay  Melnik  and  his  coworkers(S)  for 
viscous  flow  over  transonic  airfoils. 
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in  Section  2,  various  formulations  and  analyses 
are  provided  wnich  reoresent,  in  part,  a  critical 
assessment  of  different  aspects  of  the  jet  flap 
formulation  given  in  Ref.  1  ana  its  extension  to 
US3  for  incompressible  and  compressible  flows. 
Section  2.1  gives  a  full  potential  formulation  for 
the  USB  problem  which  can  be  utilized  for  future 
computational  solutions  and  as  a  framework  for  the 
more  approximate  thin  airfoil  theories  wnicn  are 
presented  in  subsequent  sections.  Candidate 
numerical  schemes  are  discussed  but  their  imple¬ 
mentation  for  the  full  potential  formulation  will 
be  reported  elsewhere.  As  a  basis  for  the  thin 
airfoil  and  small  disturbance  models,  the  thin  jet 
approximation  is  described  in  terms  of  a  systematic 
asymptotic  expansion  procedure  in  Section  2.2  for 
incompressible  flows.  The  generalization  of  these 
developments  to  compressible  *1cw  is  straigntfor- 
ward  and  therefore  not  provided.  Section  2.3  pro¬ 
vides  a  discussion  of  the  trailing  edge  -egion 
f’-om  the  viewooint  of  a  nonuniformity  of  the  thin 
jet  theory  as  well  as  the  generalized  <utta  condi¬ 
tion  for  USB.  Inherent  in  this  aspect  is  the 
geometry  of  the  dividing  streamline  at  the  trailing 
edge  which  is  a  necessary  condition  for  the  deter¬ 
mination  of  the  jet  sheet  free  boundary.  Arguments 
are  provided  to  suostantiate  tangency  to  the  upper 
surface  providing  the  jet  stagnation  pressure  is 
greater  than  that  of  the  external  flow.  To  relate 
the  internal  jet  structure  to  the  external  flow 
with  emphasis  on  boundary  conditions.  Section  2.4 
illustrates  how  asymptotic  matching  principles 
provide  a  systematic  approximation  to  the  Poundary 
conditions  for  a  problem  for  the  external  flow. 

Also  indicated  is  how  the  small  disturbance  jet 
flao  formulation  of  Ref.  1  is  .modified  with  USB. 
finally,  Section  3  gives  'esults  from  a  computa¬ 
tional  solution  based  on  the  small  disturbance 
formulation  of  the  earlier  sections,  in  this  sec¬ 
tion,  transonic  USB  airfoils  are  analyzed  and 
comoarisons  are  made  with  experiment.  ractors 
associated  with  the  discrepancies  are  considered 
and  refinements  are  oroposed  to  improve  the 
reel  ism  of  the  model . 

2.0  Emulations  and  Analyses 
Z.'  -uil  Potential  Theory 

In  rig.  1,  a  USB  configuration  is  shown. 

Before  dealing  with  small  disturbance  theory,  we 
indicate  the  aoorooriate  formulation  in  a  full 
potential  framework. 

One  procedure  to  treat  the  problem  is  to  'intro¬ 
duce  separate  potentials  *or  the  internal  'low  in 
the  jet  region  'denoted  as  Rj)  and  the  external 
-eg ion  However,  it  is  to  be  noted  that  f rota¬ 
tional  *'fow  will  result  only  if  the  entrooy  is  con¬ 
stant  across  the  jet  exit  plane  AE.  thus,  the 
mocel  will  -epresent  a  subclass  of  other  practical 
cases  wnere  vorticity  is  introduced  in  the  internal 
ducting  upstream  pf  the  jet  exit.  Hithin  the 
-estrictionj ,  the  aoorooriate  equations  are 
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where  bo  ana  b<  are  the  potentials  in  tne  regions 
Re  and  $j,  respectively.  Under  the  assumption  of 
i  rotational  i  ty ,  the  following  energy  invariants 

exist  in  R_  ana  R.: 
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const. 


(2a) 
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♦  -j-  *  const. 
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wnere  the  subscript  0  signifies  stagnation 
conditions . 


On  tne  slip  lines  AB  and  CO,  two  sets  of  bound¬ 
ary  conditions  are  recuired  since  these  lines  are 
->ee  in  the  sense  that  they  interact  with  the 
solution,  "he  #1-st  of  these  is  the  tangency  con¬ 
dition  in  which  these  lines  are  streamlines  of  the 
internal  and  external  flow,  "hus,  on  A8,  which  in 
the  coordinate  system  indicated  in  Fig.  1  is  given 
by 

y  ’  sj(x) 

the  conditions  are 


wnere  the  arguments  of  the  other  members  are  the 
same  as  the  numerator  of  the  first.  In  addition  to 
(3),  a  jumo  condition  obtained  from  Egs.  (2a)  and 
(2b)  apolies  across  AB  and  ECO.  This  -elation  is 
synonymous  with  the  fact  the  static  pressure  is 
continuous  across  each  slip  line,  so  that  across 
these  boundaries 
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holds,  where 


(*a) 


A 


e 


p(Y-l)/v 

_!a _ 


g (y-1 )/Y 


(4b) 


4c) 


Similar  conditions  hold  on  BCD.  In  sedition  to 
(3)  and  (4),  suitable  conditions  on  the  slooe  of 
AB  and  ECO  at  the  trailing  edge  A  and  C  are 
-equired  to  comolete  this  cart  of  tne  formulation, 
'hese  would  be  aetemined  by  continuity  pf  pressure 
and  flow  angle  at  these  ’ocations, sresumaoly  by 
some  iterative  procedure. 
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On  the  surface  of  the  airfoil,  a  condition  of 
the  type  (3)  applies,  i.e.,  if  the  upper  surface 
is  given  by 

y  *  fu(*) 
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A  similar  relation  holds  for  the  lower  air-'oil  sur¬ 
face.  For  the  far  field, in  R  ,  a  relation  of  the 
type 


*  *  '-jz  .  as  r2  «  x!  ♦  i1  ■*  * 

(Y  *  scaled  y) 

apolies  approximately,  where  7  is  the  circulation 
on  a  circuit  at  r  ■  »  sna  3  »  tan'2  Y/x.  Note  that 
the  circuit  includes  contributions  along  the  jet. 

The  'ar  field  in  R<  can  be  ootained  by  methods 
similar  to  that  developed  by  '•almuth  and  Murphy  in 
Ref.  5.  To  complete  the  formulation,  initial  con¬ 
ditions  at  the  jet  exit  and  a  suitaole  far  field 
downstream  are  required.  If  the  former  is  subsonic, 
Neumann  or  Oirichlet  data  are  appropriate.  If  it 
is  suoersonic,  then  Cauchy  data  are  necessary, 
downstream, data  are  used  only  if  the  flow  is 
locally  subsonic. 

At  least  three  alternatives  appear  possible  for 
the  solution  of  the  aforementioned  boundary  value 
problem.  In  t.ne  'irst,  no  maooings  are  applied  and 
the  problems  are  to  be  solved  in  tne  physical 
(Cartesian)  plane.  This  approach  has  the  advantage 
of  allowing  highly  develooed  techniques  to  be 
applied  to  the  Cartesian  fom  of  the  full  potential 
equations  (FPE).  Successful  implementation  would, 
however,  require  means  of  dealing  with  the  movable 
free  boundaries.  In  particular,  stable  extrapola¬ 
tive  and  interpolative  procedures  would  have  to  be 
develooed  to  handle  large  free  boundary  shifts 
between  grid  points,  with  minimum  error  propagation. 
In  the  Cartesian  equation  acoroach,  it  is  envisioned 
that  the  problem  would  be  solved  separately  in 
domains  in  Ri  and  R»  on  alternate  cycles  using 
successive  line  overrelaxation  ( SLOR ) .  Conditions 
in  the  other  domain  would  be  updated  using  "latest" 
data  crom  either  So.  (3)  or  (4). 

The  second  metnoo  would  use  mapping  procedures 
such  as  these  due  to  'homoson.  Another  possibility 
fs  to  ooen_uo  the  cut  physical  plane  with  the  mao- 
oing  w  »  ,7,  where  z  *  x  *  iy  and  w  »  u  *  iv.  In 
this  mapping,  the  jet  boundaries  would  project  as 
lines  /  *  Siu).  A  suosecuent  shearing 
transformation 


developed  to  effectively  capture  shocks  and  treat 
the  mixed  field  in  an  accurate  ana  computational ly 
efficient  way.  Issues  of  stability,  diagonal  oom- 
inance  would  have  to  be  resolved.  Similar  ques¬ 
tions  would  also  arise  in  connection  with  a  third 
formulation  Involving  the  stream  function  u  as  a 
dependent  variable.  In  contrast  to  the  inaepencent 
variable  case  wnere  the  free  boundary  position  is 
known  with  unxnown  ooundary  conditions,  this  would 
generate  an  unknown  boundary  carrying  xnown  cats. 

In  spite  of  the  attendant  problem  areas  in  tne 
full  potential  formulation,  it  -epresents  a  desir¬ 
able  longer  range  objective  to  treat  thick  configu¬ 
rations  ana  assess  the  accuracy  of  small  disturb¬ 
ance  'JS3  solutions.  In  the  short  term,  the  small 
disturbance  solution  nas  been  implemented  because 
of  its  simplicity  of  extension  to  three  dimensions 
and  its  ease  in  integration  with  other  inverse 
procedures  that  we  nave  developed  to  remove  sneexs 
on  unblown  airfoils.  Within  this  framework,  the 
free  boundaries  can  be  more  easily  treated  due  to 
the  appropriateness  of  transfer  of  boundary 
conditions. 

2.2  Thin  Jet  ~heory 

As  an  essential  ingredient  of  a  small  disturb¬ 
ance  formulation,  tne  iet  structure  is  developed  in 
this  section  for  purposes  of  specification  of  the 
boundary  conditions.  In  particular,  it  will  be 
shown  how  the  Spence  theory  of  Ref.  1  can  be 
derived  *>om  a  systematic  approximation  procedure. 

Referring  to  rig.  2,  a  section  of  the  jet  is 
shown.  A  curvilinear  coordinate  system  is  embedded 
in  the  jet  as  indicated.  The  lines  n  *  constant 
are  .parallel  to  a  reference  line  (the  5  axis)  which 
only  unaer  special  circumstances  ceincides  with  the 
centerline  of  the  jet.  Otherwise,  the  axis  is  a 
reference  line  wnich  is  the  centerline  of  an 
aporoximate  parallel  -low  to  be  discussed  subse¬ 
quently.  In  this  coordinate  system,  tne  lines 
5  *  constant  are  normals  to  2  axis.  In  wnat  fol¬ 
lows,  the  incompressible  case  will  pe  discussed. 

The  generalizations  to  comoressible  flow  are 
straightforward. 

Within  tne  indicated  coordinate  system,  tne 
exact  equations  of  motion  are 

Conti nui tv 

~  05  +  -rx  h  q_  *  3  (ca) 


Momentum 


.  3r<s  oh  _  -1  3D 
3n  h  gn  qn  55 


v  -  S  ( u ) 


-i  -  Momentum 


would  bring  the  '-ee  bouneary  on  the  -eal  axis. 
yet  another  naooing  procedure  would  be 
(x,y)  -  ix.ii),  where  j  ’s  the  stream  'unction  sat¬ 
isfying  the  exact  continuity  equations,  "he  advan¬ 
tage  of  these  maooings  is  that  they  woulo  "freeze" 
the  -Vee  boundar'es.  A  disadvantage  is  that  they 
would  cpmoromise  the  equation  in  the  sense  that  a 
new  numerical  methoa  would  probably  have  to  be 


dr  3Q„  3q_  di  . 
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;h  ,  j  gp 


where  h,  the  metric  coefficient  when  related 
differential  arc  length  in  Cartesian  (x,y) 
coordinates  is 

dx:  +  dy2  *  h;(5,n)d52  *  dr, 2 
h(5.n)  *  1  -  -/R(5) 
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with  9(2)  peing  the  radius  3 f  curvature  which  is 
shown  positive  in  Fig.  2. 

To  obtain  an  approximate  incompressible  set  of 
equations  prototypic  of  the  compressible  case,  the 
thin  jet  limit  is  considered.  The  characteristic 
jet  thickness  is  shown  in  Fig.  2,  wnere  the  jet 
boundary  is  denoted  as  n  *  tb(2). 

We  now  define  a  thin  jet  limit 

T*0,  2,n*  *  n/t  fixed  (7) 

where  the  boundary  layer  coordinate  m*  is  intro¬ 
duced  to  keep  the  jet  slip  lines  in  view  in  the 
limit  process.  In  (7),  the  appropriate  representa¬ 
tions  to  yield  a  nontrivial  structure  are 

qrU.n;T)  , 

-L—j - *  —  UgU.n*)  +  't-  i^U.n*)  *  •••  (3a) 


q 

TT  * 


f  •  •• 


(8b) 


where  symmetry  has  been  assumed  to  leading  order, 
and  the  jet  'low  has  been  divided  into  two  portions 
2  <  and  2  >  2t£.  with  2te  representing  the 

trailing  edge  position  and  trie  'act  that  the  trail¬ 
ing  edge  is  a  streamline  (no  jet  width  perturba¬ 
tion)  is  introduced  by  the  Heaviside  'unction  H, 
defined  as 

H(x)  *  0  ,  x  <  0 
=  1  ,  x  >  0  . 

3ased  on  the  foregoing  discussion,  the  condi¬ 
tion  that  jet  boundaries  S  are  streamlines  is 

q  •  7S  »  0 

where  q  *  (qr,qn).  Substitution  of  the  expansions 
(8)  into  thif  relation  gives: 


(11a) 

Vg('.-bg)  *  -0^(2)  UgU.-bg)  l 

(Hb) 

—  *  pn  +  tp,  +  •••  (8c) 

oUz  J  1 


where  LI  is  some  typical  freestream  velocity.  The 
orders  were  selected  to  give  the  "richest"  possible 
sec  of  equations  and,  consistent  with  this,  produce 
forcing  terms  in  the  equations  'or  the  second  order 
quantities.  These  orders  are  consistent  with  the 
massless  momentum  source  model  of  Spence. (1) 


Substitution  of  (3)  in  (5)  and  equating  like 
orders  gives  the  'ollowing  equations  'or  the 
aporoximate  quantities : 


3 


j 


^ 0 
3  38 


0 


J0  3o0 
TTTT  *  '  0^ 


3n* 


1  J_ 
T  3 r* 


(l*Vg) 


j  (9a) 

■  Ob) 
)  (9c) 
k  (19a) 


3u. 

'^2  u 

w  '’U3 

'o  F  ‘ 

317 

“0V0  n* 

u  3u0 

.  !_2 

“T-  '  T 

Jo  TT 

(10b) 


30 

3-* 


■J  (10c) 


~he  appropriate  boundary  conditions  ’Involve 
statements  -egardiog  the  'act  that  the  jet  bound¬ 
aries  are  streamlines  and  that  the  static  pressure 
is  continuous  across  the  slip  lines.  The  upper  and 
lower  slio  lines  3  and  3,  *  0  are  given  by 


V^S.bg)  *  b-j(i)  UgU.bg)  (lie) 

V-,U,-bg)  *  -b-jU)  UgU.bg)  *  SgOg 


+  b6 


(lid) 


the  other  conditions  involving  continuity  of 
pressure  are: 


PgfS.bg; 


o(u) Irl 

ext'"' 

cUJ 


2 

(12a) 

3b. 

fT  U.bg) 

(12b) 

>  present  argument 

that  the 

'ul 

external  pressure  field  is  prescribed.  In 

addition,  t  <<  £  has  been  also  implicitly  assumed. 
For  r  *  0(£),  (12b)  would  be  presumably  changed  to 
reflect  a  second  order  perturbation  of  the  outer 
flow.  In  actuality,  it  is  determined  ’Iteratively 
and  interacts  with  the  internal  jet  'low. 


Solutions 

First  Order  Theory 

We  introduce  the  cero*'*1  order  stream  function 
given  by 

3'fc  3'd 

tt  *  •  /.  ,  *  u. 

and  CTOloy  the  'allowing  transformation  r  the 

independent  va-’ables 

(?,n*l  -  U.w)  .  (13) 

Jnder  (13),  the  dif'srential  ooe^ators  mao  as 
follows: 


s^r-Tbgt-::  --^1(2)  *o  .  (^(Srg^-o) 

Sl*-*Tb0(5)*TJb1(5)H(;-2TE)  .  0 
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and  £ds .  (9)  oecome 
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"or  simplicity  and  without  excessive  lass  at' 
generality,  we  consider  a  constant  velocity  jet 
exit  (E»Q)  initial  profile  parallel  to  the  wall, 
i.e. , 


o.(0, n*) 


C  -  'jg(0,n*)  *  C 


q  (O.n*)  *  0  *  Vg(0,n»)  .  (15b) 

Equations  (14)  have  the  following  solutions  on 
aoolication  of  (15),  (11a),  (lib)  and  (12a): 


*  0 


Bq(E ,n*l  *  (1-n*)  j-  *  qfJ(S) 


Discussion 

Equations  (15)  describe  a  parallel  ‘low  jet. 

The  total  jump  in  pressure  across  the  jet  ‘-om 

(16c)  is 

Cp0J  *  q(E.l)  -  o(5. *1 )  *  2C*/R  (17) 

which  agrees  with  the  Soence  model.  It  should  be 
noted  that  in  contrast  to  the  latter,  no  assumption 
regarding  irrotational ity  is  required  to  obtain 
(17),  in  contrast  to  the  results  of  previous 
workers.  The  radius  of  curvature  of  the  jet  is 
approximately  9  upstream  of  the  trailing  edge  which 
in  turn  is  approximately  given  by  that  of  the  blown 
upper  surface.  Downstream  of  the  trailing  edge, 

R  is  determined  ‘rom  applying  (17)  to  the  determi¬ 
nation  of  the  flow  outside  of  the  jet.  Upstream  of 
the  trailing  edge,  the  wall  pressure  is  determined 
oy  (17)  since  R  is  known  and  is  given  by 

P(E>-1)  *  *  9,j(E)  •  (18) 

Eecond  Order  Theory 

Equations  (10)  have  been  solved  subject  to  the 
boundary  conditions  (He),  (lid),  (12b)  and  the 
initial  condition  at  the  jet  exit  -;*0 

u,(0,-.»)  •  0 


to  give  the  following  solution,  wnere  the  first 
order  solution  ;16)  is  assumed: 

(  ,  q'(r) 

i  VV1 

!  *  +  ,  0  <  Z  <  Ete  (19 

I  #  q 1 

'  V1  *  V1  "(T  E1*'!*-  *  uQ(  1/R)  * 
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-  [z—'-rj  ■  s> 

vl(hTr*,-’ )  *  -ql_(0+)/Ug 


’  ^  ’TE 


■^(^g-.-l)  *  0  (19d) 

u1  *  u0(l-n*)[R(0)-‘  -  R‘!]  ♦  ( q.j ( 0 )  -  q.J/Ug  (19e) 
b,  -  2Cqu  -  qu(0)2/u§  (19f) 

.  2u5  f n*‘  ,  ll  ,J0  r,  ,  3  „  ll 

.T‘"  IJr^T  [_2n  'I'  ?J 

-  ¥■  Ci./°)  -  qu-  ■  os?) 

Besides  u g  -  0,  another  nonun iformi ty  is 
evident  in  (l?a)  and  (19c)  and  is  associated  with 
the  factor  q,'(;j  as  E  -  E-£,  if  the  latter 

approaches  ».  As  we  will  see,  this  occurs  ‘or 
trailing  edge  angles  less  than  SO0. 

2.3  Trailing  Edge  3enavior 

Incompressible  flows 

Defining  the  complex  potential  ‘unction  as 
F( 2)  *  a(x,y)  +  i'p(x.y) 

where  (x,y)  is  the  local  coordinate  s>j*am  shown 
in  Fig.  3  and  z  »  x+iy,  the  local  "corner  ‘"ow' 
solution  to  within  a  dimensional  multiplicative 
constant  in  the  lower  external  region  A'DB  with 
i  as  the  dividing  streamline  angie  shown  is 
given  by 

F  .  8-M-r+v)/a  ,r/S  (20) 

which  implies  that  the  square  sf  the  resultant 
velocity  q  is 

| F‘ ; 2  «  q2  "«  r2i?  *  ')  .  (20‘) 

To  obtain  the  necessary  initial  conditions  ‘or 
determination  of  the  jet,  the  implications  of  con¬ 
tinuous  pressure  and  flow  angle  across  the  dividing 
streamline  near  the  trailing  edge  were  studied  ‘or 
incompressible  ‘low.  Referring  to  Fig.  3,  with  the 
dividing  streamline  denoted  as  OB,  we  signify  the 
trailing  edge  angle  as  v,  the  angle  that  OS  maxes 
with  the  upper  surface  AO  as  -  and  that  with  the 
lower  surface  as  3.  In  what  ‘ollows,  we  identify 
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each  of  nine  possible  cases  Chat  can  occur  with 
the  -'ol  lowing  shorthand  notation  exemplified  as 
f o 11 ows : 

LG  •  j  <  T,  3  >  T 
£5*^*1,  3  >  «r 
(3  ■  2t  -  j  -  v)  . 

The  'truth  taole"  indicating  the  nine  cases  is 

r  LS  L£\ 

(  GL  GG  GE 


?  *  P0  -  ~y~  *  p  *  ?0 

where  *  is  the  upper  stream  soeed  at  0.  To 
satisfy  (21),  the  "slip"  velocity  S'  is  thus 


n+  •  V- 


2(pn  -  Pn) 


Remarks 

If  pt  <  pi,  then  IE  is  the  only  possibility 
with  u  J 


Introducing  the  notation  in  which  *  subscripts 
and  superscripts  signify  conditions  above  as  and 
-  subscripts  and  superscripts  denote  points  below 
this  line,  we  now  examine  each  of  the  possibili¬ 
ties.  In  wnat  follows,  let  pq  be  tne  stagnation 
pressure  and  o  be  its  static  value.  Except  where 
otherwise  noted,  we  will  also  assume  that 

O0  >  ■ 

All  Possibilities  Involving  G  Except  GG 

These  must  be  ruled  out  on  the  basis  that  on 
the  G  side,  p  *  wnicn  is  unphysical  on  its  own 
merits,  furthermore,  p  would  equal  pg,  the  stagna¬ 
tion  value  on  tne  L  side,  and  equality  of  pressure 
across  tne  slip  line  would  therefore  "be  impossible. 


For  a  hypothetical  case  in  which  in  a  real 

physical  flow  (fluctuating  above  and  below 
equality),  then  a  tri-stable  configuration  could 
evolve  which  would  oscillate  between  LL,  l£  and  EL. 

Monun if oral  ties  of  Second  Order  Approximation 

From  (20),  we  find  that  the  other  nonunifpraity 
referred  to  following  Eq.  (19)  is  associated  with 
the  following  behavior 

qu  ^  r^2,/8^’3  as  r  -  0 

where,  we  distinguish  the  -allowing  possibilities 
as  i  -  iT£: 


This  must  be  excluded  since  it  would  imply 
3  *  »  >  2-  which  is  clearly  geometrically 
impossible. 


This  implies  that 
*  ♦ 


To  realize  pressure  equality  across  the  slip  lines, 
this  must  be  ruled  out  as  impossible  for  all  cases 

except  when  o^  *  pi.  This  also  is  a  corollary 

of  the  Kutta  condition  applicable  to  unolown 
configurations. 


(1) 

qu  *  *  • 

3  >  or 

T 

-  <  T 

(11) 

q;*o  . 

3  <  Ft  or 

u  >  j 

(ill) 

q,j  finite 

*0,3* 

T  or  'J  *  j 

Case  (i)  is  the  most  practical  situation  and  will 
necessitate  an  inner  solution  ‘or  the  transition 
layer  to  join  the  wall  and  ‘ree  jet  flows.  9efore 
considering  this,  we  briefly  investigate  the 
vorticity  wnicn  can  be  snown  to  be  given  by 

M  *  t  ojq  +  T*  uij  *  •••  ,  wnere 


'his  is  associated  with  v  *  G,  i.e.,  a  cusped 
trailing  edge  or  flat  plate  for  *  Pg. 


Jere,  p"  <  Pg  <  pg.  However,  p*  *  Pg  in  this 
case  wnicn  makes  aoui 1 ibration  of  the  static  pres¬ 
sures  o"  and  p  ‘moossible  across  the  slip  line  A8. 


+  "his  is  the  only  /‘able  possibility.  Here, 
P  <  Pa.  In  ‘act,  by  Sernoullf  and  slip  line 
pressure  continuity,  we  have 


‘U1  ug  C 

*  V?  *  T  ’  •  ’  C3nstant  *  °  ' 

In  soite  of  the  constant  initial  ‘low,  a  non-cero 
vorticity  is  thus  introduced  by  the  body  curvature. 
3y  contrast,  a  ootential  vortex  over  a  circular 
cylindrical  sur-'ace  would  have  had  a  linear 
initial  profile  to  produce  an  irrotational  flow. 

Inner  or0Q)CT 

Anticipating  large  gradients  near  ‘  *  ;T.  we 
consider  the  following  asymptotic  expansions" /al id 
in  the  inner  limit 

v*  *  (;-'rr)A(f )  ,  “*  *  *7-  fixed  as  T - 0  23) 
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where  C  has  seen  assumed  do  be  unity  in  (15): 


where  F  is  an  arbitrary  Function  or  integration. 
The  system  (28)  is  equivalent  to  the  following 
scalar  equations 


—  *  s(t)u *(C*,n*)  *  ••• 


—  *  3(x)v*(i*,n*) 


^  *  P*(i*,n*)  +  ••• 

ell2 

where,  p  is  a  reference  pressure  and 


h  *  1  - 


Wl *7 


Mere,  R*  represents  the  centerline  radius  of 
curvature. 

To  obtain  the  richest  continuity  equation 
assuming  R*  >  Q(t),  we  out 


to  give 


1H1  +  2*1 , 0 

In*  * 


Roughly,  the  matching  condition  for  j*  is 


-  0  t 


u*  -  0  . 


For  v*,  we  introduce  the  intermediate  variable 
v;  ■ 

Matching  the  inner  and  outer  solutions  by  writ¬ 
ing  each  in  intermediate  variables,  we  assert  that 

3(TMr  vn*)  -  ^  bo(vVuow  ” 

where  the  double  arrow  signifies  asymptotic  equal¬ 
ity.  Equation  (25)  implies  that 

bj(01  »  v*(~»,r*)  *  constant  *  0 


3  '  /t  ,  s/1  ‘  1//F  . 

Substituting  the  expansions  (24)  into  the 
momentum  equations,  the  remaining  equations  in 
the  distinguished  limit 

\*T,3*/t,a*/t  ( 


AP*  =  21=1  *  *  0  (29a) 

SC*2  3n’2 

Av*  *  {-  1 /R*) 1  .  (29b) 

The  appropriate  boundary  value  proolem  for  (29) 
involves  specializations  at  boundary  conditions 
already  discussed  and  conditions  obtained  *rom 
matching  at  the  upstream  and  downstream  boundaries 
of  a  rectangular  domain.  The  jet  boundaries  are 
not  free  to  this  order,  oeing  at  n*  *  ;1 .  A  com¬ 
plete  solution  of  the  problem  subject  to  these 
conditions  is  in  progress.  As  an  interim  step,  a 
special  form  of  an  "inner  inner"  solution  is  indi¬ 
cated  *or  the  jet  flow  immediately  aoove  the  trail¬ 
ing  edge.  To  illustrate  the  singular  benavior, 
vorticity  associated  with  a  source  term  in  the 
Poisson  equation  for  u  obtained  from  (28)  and  (25) 
is  suppressed  by  a  choice  of  suitable  jet  exit 
conditions.  In  addition,  tne  ’•ignt  hand  side  of 
(29b)  is  assumed  to  vanish  by  1 inearization  about 
;  =  iTE  in  The  subscale  assumed. 

For  purposes  of  the  analysis,  the  configuration 
in  Fig.  3  is  considered.  Here,  ^  =  n,  representing 
the  EL  arrangement  which  was  proven  previously  to 
be  appropriate. 

On  applying  Bernoulli’s  theorem  on  the  uooer 
and  lower  sides  of  the  trailing  edge,  and  employing 
the  condition  that  the  pressure  is  continuous 
across  the  dividing  streamline  38  we  obtain 
using  (20’1 


if  e  u(x,Q+)  * 


30  *  p0 


/(H 


as  x  -  0 


2S2  Pn  -  Pr, 


where  o  is  the  density. 

The  solution  for  the  jet  flow  in  tne  immediate 
neighborhood  of  the  trailing  edge  satisfies  the 
approoriate  boundary  value  problem  for  the  upper 
half  plane  y  >  0  (Fig.  4)  in  wnich  flow  angle  and 
pressure  are  matched.  Here,  continuity  of  flow 
turning  signifies  that  in  some  "inner"  ’•eg ion 

v(x,3)  »  3  ,  x  <  0  (31a) 

v(x,0)  *  0  ,  x  -  0  (31b) 


-V*  1  ^  l  ;Q« 

3£*  *  R’  °  ”  c  in* 
3u*  .  1  30* 


Equation  7 2Sb )  can  be  integrated  directly  to  give 
the  linearized  Bernoulli  relation 


u(x,C+)  =■  A  *  3xa  (32) 

where  A  and  3  are  constants  determined  from  (31). 
The  appropriate  harmonic  functions  satisfying 
(31)  and  (32)  are 

F’  (z)  s  a-  *  a,e''TCl  z*  (33a 


u*  «  .  •  F(r*) 
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u  *  a^a^r1  cos  x( 3-m)  *  Re  F’(z) 


(33b) 


for  which 


Mith  '.he  usual  isertrooic  relation 


u(x.O) 

*  a^  *  a,x3*  ••  • 

.  x  *0 

(34a) 

u(x,'J) 

*  a„  *  a,x*  css  am  *  •  •• 

U  i 

1  X  r  0 

(34b) 

v(x,3) 

*  3 

.  <  *  0 

(34c) 

v(x,0) 

*  *v:‘3in  3" 

,  X  r  0  . 

(340) 

'he  shape  of  the  sliD  line  is  determined  by  apoly- 

i ng  Eos.  (34)  in  c.ne  : angency  relation,  denoting 

the  jet  ooundary  by 

J  *  y  -  3(x)  *  0 

then  the  flow  tangency  relation  (u,v)  •  73  *  0 
implies  that 


3’<x)  *  F'!y.3(x) 


-1 

a^  *  a^3  cos  am  a0 


as  x - 0  (35) 

for  3  <  t  . 


Hook  of  dividing  Streamline  -  Discussion 
Equation  (35)  implies  that 

3" (x)  *  -a  x3*1  . 


3"(x)~»  as  x  -  0  if  v  <  y 

in  spite  of  the  fact  that 

3'  (x)  -0  as  x-0  , 

i.e.,  a  'hook"  of  infinite  curvature  but  with  zero 
slooe  develops  at  the  intersection  of  the  dividing 
streamline  and  the  trailing  edge.  This  fact  has 
significant  conseauences  'egarding  the  generation 
of  lift  of  blown  airfoils.  In  this  context,  the 
classical  Soence  solution  'ails  to  treat  this 
detail  since  it  is  cast  in  a  thin  air*oil  theory 
framework.  It  gives  a  logarithmic  singularity  at 
the  trailing  edge  in  contrast  to  that  of  the  solu¬ 
tion  above,  'he  '•soresentation  of  3  near  x»0  from 
(35)  provides  the  initial  conditions  for  the  deter¬ 
mination  of  the  lower  boundary  of  the  jet  and  the 
vorticity  it  carries,  "his  in  turn  determines  the 
running  circulation  and  the  total  lift.  A  systema¬ 
tic  matching  procedure  to  Join  this  solution  to  the 
outer  -low  and  achieve  a  consistent  numerical 
'ormu’ation  nas  aooarently  not  been  -•eoorted  in  the 
literature,  even  'or  the  unblown  incompressible 
case  cor-esoonding  to  thin  airfoil  theory,  i.e., 
epr  x<  1 . 

Icmoress''ble  'raii'rq  Edces 

Consider  again  the  configuration  of  cig.  3. 

Mere  we  analyze  first  the  case  where 


30  ’  ao  • 


the  local  Mac.n  numoer,  M ,  can  adjust  in  a  continu¬ 
ous  way  in  order  that  the  'low  decompresses 
smoothly  from  3  to  0.  At  0,  M  is  single  valued**^ 
and  adjusts  itself  such  t.nat  the  pressure  p* 
eouals  ot,  in  accord  with  the  <sentrcpic  relation, 
assuming Jthat  the  flow  stagnates  at  13  on  the  lower 
sice.  The  only  way  this  can  be  realized  is  with 
an  EL  configuration.  3y  contrast,  a  oiscontinuous 
transition  can  occur  and  is  illustrated  schematic¬ 
ally  in  Fig.  5a,  and  is  associated  with  a  GL 
arrangement.  In  some  respects  tne  configuration 
resembles  flow  over  a  solid  wall  expansion  corner 
consisting  of  expansion  *an  interacting  with  a 
sonic  line  *rom  the  comer.  If  the  solid  configu¬ 
ration  were  representative  of  this  flow  with  the 
free  boundary  slip  line,  compression  waves  would 
reflect  of*  the  sonic  line  anq  -arm  a  snocx 
envelope  which  would  oe  necessary  to  decompress 
the  flow  from  an  overe.xoanaeo  value  oelow  the  crit¬ 
ical  p+  to  the  oj  level.  Additional  reflections 
can  be  produced  from  tne  upper  slip  line  30. 

This  discontinuous  transition  ’eads  to  a  multival¬ 
ued  pressure  at  0.  rhe  continuous  and  discontinu¬ 
ous  processes  are  illustrated  schematically  in 
Figs.  5b  and  5c,  respectively.  In  F<g.  5b,  the 
-ecompression  takes  alace  on  the  line  30.  In 
Fig.  5c  it  occurs  on  00'C”0‘"  signifying  the 
confluence  of  multiple  states  at  3.  Mere,  the 
dashed  line  element  0"0"'  signifies  a  shock  jump. 
Experimental  data  strongly  suggests  that  as  in  the 
incompressible  case,  the  EL  configuration  is  the 
most  probable  situation, since  on  a  qualitative 
basis.it  represents  a  path  of  least  -esistance. 
Presumably,  a  more  r-gorous  argument  to  suooort 
this  conjecture  would  -ely  on  some  sort  of  stabil¬ 
ity  analysis. 

For  pi  =  ol,  the  dividing  streami'ne  would 
again  bisect  the  trailing  edge  since  the  '1ow  in 
the  immediate  vicinity  would  be  incompressible  ana 
the  reasoning  in  the  previous  section  would  apply. 
For  the  improbable  case  of  pt  <  a;  the  uE  configu¬ 
ration  would  be  applicable  as  ‘‘orJthe  incompressi¬ 
ble  situation. 

3.4  incompressible  Small  Deflection  '-eo rv 

3ecause  of  its  potential  value  eor  simpi  i-'-'eo 
treatment  of  the  trailing  edge  neighborhood  ana 
matching  with  the  outer  flow,  we  consider  ;n  this 
section  a  small  deflection  specialization  of  the 
previous  thin  jet  theory,  -or  surooses  of  illus¬ 
trating  the  theory,  we  -evert  back  to  incompressi¬ 
ble  flow.  >ie  consider  a  section  of  the  jet 
far  away  ‘>om  the  trailing  edge  -egion  depicted 
in  Fig.  5.  "he  jooer  and  lower  slip  lines  a-e 
given  sy 

J  *  y  -  it>u(x)  *  0  (36a) 

J,  *  y  -  5bz(x)  *  0  (36b) 

where  the  slooe  of  the  slip  line  is  assumed  to  be 
of  the  same  order  of  the  characteristic  thickness 
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ratio  o  or  incidence  a  of  one  airfoil.  Moreover, 
(36)  imolies  that  the  jet  thicxness  is  also  of  the 
same  order. 

We  assume  asymptotic  expansions  for  the  veloc¬ 
ity  vector  o  and  pressure  ?  can  be  written  in  the 
form 


4^-*  [y  SjUgU.yJ  +  ^u 

+  [3o(5)vo+"-]j' 

P 

T"  *  <0(5)P0  +  <1(i)P1  +  --- 


J  -  3b. (x) 

for  x  *  x,y  *  - A -  -ixed  as  i-0,  where  'J  is 

the  freestream  velocity,  :  is  the  density,  and  P^ 
is  the  ambient  pressure.  The  condition  tnat  the" 
lower  slip  line  is  a  streamline  of  the  jet  flow 
reads 

0  •  ~J.  *  0 

wnich  to  dominant  order  implies 

v„(x,0)  *  -j-2  b  1  ( x }  (38a 

u  -0  <■ 

•which  provides  a  non-trivial  case  only  if 


Case  (ii)  appears  to  be  the  most  interesting  of 
the  three  since  A)  leads  to  constant  vertical  and 
horizontal  velocity  along  streamlines  which 
restricts  the  free  surface  snaoe  and  (iii)  leaos 
to  no  pressure  difference  across  the  jet  to  this 
order  which  is  of  limited  practical  interest,  -or 
case  (ii),  we  introduce  tne  stream  function 
mapping 

(x.y)  -  (x,p)  (41 ) 

in  which 


'O'  ,  v.  *  -  •!) 
y  0  x 


lx  *  U 


3y  J0  3 'h 


3asad  on  (38b),  a  similar  relation  is  obtained 
on  the  uoper  slip  line  which  is 

VgU.bj-b^)  *  bj(x)  .  (38c 


3  s  Uw  -4  . 

0  3x 

Under  (40)  with  (41),  Egs.  (39)  become 


x  Momentum 


On  substitution  of  (37)  and  (38b)  into  the  exact 
incomoressible  Euler  eouations,  the  following 
system  evolves: 

x  Momentum 


3u0  =  '  ~ 


Momen  turn 


52’o  °y 


Continui ty 


V  *  3 


v  Momentum 

%0m  ^0 
IT  a  3-1) 


Continuity 


te  •  *) 


V"0  £  • 3 


The  solution  of  (42)  subject  to  the  oounoary 
conditions  is 

uQ  ■  Ub) 


Vg  *  bj(x)Ub) 


v  »  3  .  k  1  3  n  s  1 1  3  ^  .. 

3x  n  S1  3?  *  3  0  3x  0  f)  ■ 

On  inspection  of  the  system  (33),  we  note  the 
following  oossibilities 


p  -  o(x)  *  -b'Ux)  f  in'  .  (43c) 

u  -  j  ^ 

Here.  o4(x)  and  b,(x)  are  the  pressure  and  slooe 
of  tne  lower  slip  line  of  the  jet.  Eos.  (43) 
indicate  that  the  Horizontal  velocity  comoonent  is 
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convectaa  along  streamlines  from  a  prescribed 
initial  distribution  at  some  upstream  station.  If 
this  distribution  is  constant,  i.e.,  if 

dg CO.*)  *  C  , 
then  (43)  specializes  to 


To  preserve  the  structure  of  the  secona  argument 
of  the  right  hand  side  of  (45),  we  assert  that 

5  -  v  . 

Furthermore,  we  obtain  the  most  general  boundary 
condition  in  (45)  if  we  let 


uQ  *  C 

Vg  *  C  b*(*) 

■u  *  y  -  b4(x) 

Pu  -  p  *  *  C  bj(x)  (ij»u  -  ’ll) 


(44a) 

(44b) 

(44c) 

(44d) 


U  *  5 2  . 

Now  if  A  *  32/3x2  *  3l/3y2,  we  have  that 
AS  *  0 

which  implies  that  one-dimensional  flow  given  oy 

Vy*  *  0  (46) 


b  -  b,  *  const. 

'J  4 

Note  that  Egs.  (44)  are  consistent  with  Eqs.  (16). 
In  the  present  small  deflection  context,  we  go 
beyond  the  analysis  associated  with  (16)  and 
'patch"  to  the  outer  flow  to  determine  the  slip 
line  deflection  functions  b,j,bj(x)  and  the  gauge 
functions  ng  and  <n.  The  term  "patch"  is  used  in 
contrast  to  "match’  in  the  sense  that  conditions 
are  satisfied  at  a  fixed  boundary  at  a  finite  dis¬ 
tance  in  the  former  (i.e.,  the  slip  line)  in  con¬ 
trast  to  the  latter  where  they  are  satisfied  in  an 
asymptotic  manner,  generally,  at  an  infinite 
distance. 

Patching  is  achieved  using  a  "blending  layer" 
at  a  vertical  distance  of  0(6)  from  the  slip  line 
and  airfoil.  The  latter  will  not  enter  into  the 
present  discussion  but  will  be  considered  in  con¬ 
nection  with  the  trailing  edge  behavior.  The 
blending  layer  will  also  be  used  to  validate  the 
usual  "aylor's  series  transfer  of  boundary  condi¬ 
tions  employed  to  define  the  outer  flow.  A  treat¬ 
ment  of  similar  blending  layers  is  discussed  in 
Cole(')  in  connection  with  incompressible  flows 
around  unblown  bodies  of  revolution.  However, 
some  different  '‘eatures  arise  in  the  present  con¬ 
text.  In  accord  with  the  previous  remarks,  we 
consider  the  shaded  '•eg ion  shown  in  Fig.  7.  We 
consider  the  entire  external  flow  “Meld  Pe  in 
Fig.  1  and  in  particular  the  blending  layer  as 
irrotational .  In  anticipation  with  matching  with 
the  outer  flow,  we  assume  that  the  blending  layer 
is  a  perturbation  on  the  freestream.  Accordingly, 
the  asymptotic  expansions  for  the  potential  i  and 
pressure  are 

g*  x  +  u($Mx.y»)  +  •••  ,  x,y*  fixed 


where 


y*  *  y/\i(i)  . 


Introducing  the  conoition  that  the  slip  line  is  a 
streamline  of  the  external  flow,  gives  on  the 
uooer  line  J  -  0  the  relation 

<i^VTV),(**bi*l)  *° 

vM  su(x))  '  £km  •  (1S 


or 


f  *  a(x)y*  +  b 

*  b^(x)y*  +  b(x)  ,  (b^  *  b’ ) 

by  (45).  A  similar  result  for  the  lower  blending 
layer  is 

V  *  bly*  b(x)  . 

i 

Now  the  well  known  outer  exoansion  for  the  region 
outside  the  "inner"  blending  layer  is 

■j  *  x  +  5d>(x ,y )  +  ••• 
x,y  fixed  as  6-0  . 

Matching  of  the  inner  and  outer  expansions  can 
oe  achieved  by  introduction  of  an  intermediate 
variable 


where  5  <<  n  «  1 .  The  matching  condition  is 
11mn  5  J3  [*outer(x-V*4inner<x’V] 


(47) 


where 


s 

outer 

~nr- 


4 . 

inner 


¥ 


x  *  S$(x,n.y  )  -  ••• 

x  +  6^i(x,0*)  *  -y,sy(x 

,0+)  +  •••]"■  i)(52) 

(4Sa) 

u. 

-y,  1 

x  *  6®(x,0+)  +  5 2 ! b ' ( x ) 

L  4 

*  b  ( X )] 

(48b) 

where  the  second  term  of  (18b)  has  been  interposed 
for  matching  and  can  be  construed  as  an  additional 
element  lumped  into  b(x)  of  the  dominant  inner 
problem  associated  with  (46).  Application  of  (47) 
and  a  similar  argument  for  the  lower  blending 
layer  gives  *or  matching  that 

i,.(x,0+)  *  b'(x)  *  j  (x,0-)  *  b:(x)  (49) 

y  u  y  v 

which  is  the  condition  anticipated  crpm  tangent 
*low.  Determination  of  b(x)  in  (48b)  depends  on 
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the  inner  limit  of  the  outer  representation  of  the 
flow  field  as  a  line  vortex  extending  to  downstream 
infinity  whose  local  strength  is  proDortional  to 
the  lump  in  pressure  across  the  jet 
[?x]  =  Jx(x.0+)  -  Sx(x.O-) .  The  subsequent  argument 
should  be  similar  in  some  respects  to  that  associ¬ 
ated  with  the  finite  line  source  discussed  in 
Ref.  7.  In  this  connection,  the  details  of  the 
determination  of  b(x)  will  be  discussed  elsewhere. 


To  determine  the  gauge  functions  <g  and  sg 
appearing  in  (37),  the  static  pressures  of  the  jet 
are  equilibrated  to  those  of  the  blending  layer  at 
the  slip  line.  From  the  interposed  term  in  ( 48b ) 
and  Bernoulli's  relation  in  the  blending  layer, 
the  external  pressure  ?a  is  given  by 


(50) 


in  both  inner  and  outer  limits,  where  on  the  slip 
1  ines 


Pe(x,bu(x))  *  -  ix(x,0+) 


(51a) 


External  F^ow 


For  distances  large  compared  to  the  jet  width, 
the  Mne  structure  of  tne  jet  is  important  only 
insofar  as  it  provides  matching  conditions  to  the 
irrotationa!  flow  field  outside  Itself.  In  incom¬ 
pressible  Mow,  this  external  "outer"  flow  can  be 
determined  by  thin  airfoil  theory.  At  transonic 
speeds,  small  disturbance  theory  is  appropriate 
for  this  region.  Details  of  the  asymptotic  match¬ 
ing  proceaure  have  been  discussed  for  incompres¬ 
sible  flow  in  the  previous  section.  Based  on  these 
developments  and  the  earlier  ones  for  arbitrary 
deflection  thin  jets  in  Section  2.2,  the  boundary 
conditions  for  the  outer  flow  in  the  incompressible 
and  transonic  cases  for  the  jet  flap  and  upper  sur¬ 
face  blowing  are  now  indicated. 

Jet  Flap 

Referring  to  Fig.  3,  the  equation  of  an  airfoil 
can  be  given  as 

y  a  jf(x) 
and  the  jet  is 


o9(x,bl(x))  *  -  ox(x,0-) 


(51b) 


y  *  5g(x) 


From  (51)  and  continuity  of  pressure  across  the 
slip  line  it  is  obvious  therefore  that 

*  5  (52a) 


and 


Ug  *  S*4  ,  30  *  S’*  .  (52b) 

Furthermore,  3ernoulli  also  implies  that 


again  in  agreement  with  tos.  (15)  for  the  jet  cen¬ 
terline  of  raoius  of  curvature  0(a"').  Note  that 
the  vorticity  uy-vx  of  the  general  solution  (43) 
is  t’tS-'a) ,  although  in  the  special  case  (44)  it  is 
This  fact  appears  net  to  affect  the  result 
(44d)  which  is  in  agreement  with  the  Spence  rela¬ 
tion  obtained  under  the  assumotion  of  irrotational 
potential  vortex  flow  for  a  jet  element. 


where  5  is  the  thickness  ratio  of  the  airfoil,  f  is 
the  upper  or  lower  surface  and  involves  the  angle 
of  attack  which  is  of  the  same  order  of  5.  Con¬ 
sidering  a  small  disturoanes  approximation,  we 
obtain 


tx 


r  s 

3/2 

r  , 

(j  *  S‘g '  * j 

W  .'o  ■'  >  X  ' 

577  '9 


R*1  * 


Letting  the  "outer"  expansion  pressure  coefficient 
be  represented  as 


P-P, 

cUi 


*  5P(x,y)  + 


then  by  virtue  of  a  generalization  of  (42) 


[P(x,0)j  *  -  Cj9"(x)  •  -  2[jx] 


(53) 


It  is  significant  to  note  that  a  oonuniformity 
occurs  near  the  trailing  edge  point  C  in  Fig.  1. 
The  order  of  the  oressure  oerturoations  in  (51b) 
are  incorrect  near  the  trailing  edge  stagnation 
region  of  the  flow  below  the  slip  line,  i.e., 


Since  oressure  is  to  be  matched  across  the  slio 
line,  this  ’’molies  a  corresoonding  change  of 
orders  in  the  jet  and  hence,  (52a)  becomes 
invalid.  This  feature  as  well  as  the  blending 
layer  structure  are  anticioated  as  important 
asoects  of  the  transonic  orcblem,  and  are  under 
investigation.  As  indicated  previously,  local 
expansions  are  required  to  deal  with  the 
nonunifomities. 


where 

Cj  s  (o  JT  q|dn^/oU’  -0(1)  (54a) 

and  i  is  a  perturbation  potential.  Note  that  (53) 
can  also  be  obtained  from  (44d). 

Equation  (53)  is  the  relation  usea  in  conjunc¬ 
tion  with  the  jet  tangency  relation 

Sy(x,0)  •  g 1 ( x )  (54b) 

and  the  airfoil  boundary  conditions  to  determine 
the  external  Mow  field.  These  relations  coincide 
with  those  derived  by  Spence,  'hey  can  be  general¬ 
ized  *or  transonic  flow  by  placing  the  c  ^nside 
the  Integrand  in  (54a). 


Upper  Surface  31owinq 

To  treat  conditions  on  the  blown  part  of  the 
airfoil  Eq.  (53)  can  be  applied  by  approximating 
the  radius  R  by  (f")*;  to  obtain  the  wall  pres¬ 
sures,  and  using  the  airfoil  ano  jet  boundary 
conditions  to  determine  the  upper  slip  line  jet 
pressures. 

rrom  the  arbitrary  deflection  thin  jet  theory 
derived  in  Section  2.2,  it  can  be  seen  that  rota¬ 
tional  flow  produces  the  same  pressure  jumps 
across  the  jet  in  the  dominant  approximation  as 
the  irrotational  Spence  models.  Correspondingly, 
it  can  be  shown  that  to  within  factors  involving 
the  density,  qualitatively  similar  results  are 
ootaineo  for  transonic  flow.  Another  important 
aspect  of  the  asymptotic  representations  derived 
here  is  that  they  lead  to  higher  approximations 
'or  the  structure  of  the  jet  and  external  flow 
whicn  can  be  systematically  obtained.  Finally, 
the  analytical  solutions  described  above  allow 
the  systematic  assessment  of  the  effects  of  initial 
vorticity  and  skewness  wnich  are  inaccessible  to 
other  theories. 

3.0  Results  and  Discussion  -or 

Transonic  Jpper  Surface  3)owing 

A  successive  line  overrelaxation  (SLOR)  scheme 
within  a  Karman  Suderley  framework  has  been  used 
to  compute  the  flow  field  over  an  upcer  surface 
blown  airfoil.  On  the  blown  portion,  the  jump 
conditions  across  the  jet  are  determined  by  the 
asymptotic  results  given  in  previous  sections, 
i.e.,  Eqs.  (53)  and  (54b).  Providing  that  the 
region  is  not  too  close  to  the  jet  exit  or  trail¬ 
ing  edge,  the  streamwise  gradients  can  be 
neglected  in  the  entropy  and  velocity  component 
parallel  to  the  wall.  Away  'rom  these  regions, 
the  pressure  gradient  perpendicular  to  the  stream¬ 
lines  is  balanced  by  centrifugal  'orce.  For  the 
region  near  the  jet  exit,  these  assumptions  become 
invalid.  Here,  the  scale  of  the  gradients  in  the 
streamwise  direction  become  important,  principally 
due  to  the  influence  of  wave  interactions  with  the 
sl'j  line.  Similar  fine  structures  occur  near  the 
trailing  edge  where  the  flow  can  stagnate  on  the 
unblown  side,  depending  on  the  ratio  of  the  stag¬ 
nation  pressure  of  the  jet  to  the  ambient  stagna¬ 
tion  value.  For  incompressible  flow,  the  previous 
sections  have  discussed  the  tri-stable  equilibrium 
at  the  trailing  edge  corresponding  to  the  value  of 
the  stagnation  pressure  ratio,  which  leads  to  the 
dividing  streamline  leaving  tangent  to  the  upper 
surface  if  this  is  greater  than  unity  ("EL"  con¬ 
figuration).  Consistent  with  the  previous  discus¬ 
sion,  the  appropriate  generalization  to  transonic 
flow  was  assumed  also  to  be  EL  for  a  single  valued 
pressure  without  a  shock  in  that  location.  *his 
assumption  has  been  altered  to  assess  the  sensitiv¬ 
ity  of  the  flow  to  the  dividing  streamline  angle, 
in  this  connection,  surface  pressures  'or  the 
dividing  streamline  bisecting  the  trailing  edge 
angle  (as  ft  would  in  incompressible  flow)  were 
compared  with  those  'or  the  EL  arrangement.  3ased 
on  these  studies,  significant  differences  are 
ant4cioaced  only  'or  large  incidences  and  trailing 
edge  angles. 

*ypical  results  obtained  'rom  the  computa¬ 
tional  model  are  shown  in  Fig.  9  in  which  the  flow 
over  a  thick  airfoil  designed  at  Rockwell’s 
Columbus  Aircraft  Division  (CAD)  was  analyzed  with 


the  SLOR  code.  Here,  the  pressures  for  various 
values  of  the  blowing  coefficient  Cj  are  comparec 
against  those  'or  the  unplown  case  at  a  'reestream 
Macn  number  Ma  =  3.733,  and  angle  of  attack 
x  *  3".  Substantial  lift  augmentation  is  evident 
for  blowing.  Also  evident  is  the  associated  rear¬ 
ward  motion  of  the  shock  with  increased  blowing 
and  sectional  loading  as  if  the  incidence  is 
increased. 

Further  parametric  studies  are  provided  in 
Fig.  13  which  indicate  the  effect  of  parallel  dis¬ 
placement  of  the  slot  xj  (in  units  of  the  chord), 
on  the  chordwise  pressures.  Three  positions  of 
the  slot  xj'0.5,  0.55  and  3.3  are  shown.  No  sys¬ 
tematic  trend  in  the  blown  pressures  is  exhibited 
on  this  air'oil  with  downstream  slot  movement  for 
fixed  Cj.  Evident  nowever  is  a  slight  intensifi¬ 
cation  of  the  terminating  snock  with  slot  down¬ 
stream  motion  although  its  position  remains  .inal- 
tered.  Despite  the  limitations  of  the  model  to 
describe  the  'ine  structure  of  the  jet  exit 
”•69100,  a  small  suction  pea*  which  has  been 
observed  in  experiments  is  exhibited  in  tnis 
vicinity  for  x,*3.5.  In  Fig.  11,  the  correspond¬ 
ing  increase  in  li't  coefficient  Ci_  with  slot 
downstream  movement  is  also  shown  as  well  as  the 
increase  in  the  size  of  the  supersonic  region. 

In  Fig.  12,  the  increase  of  lift  with  blowing 
coefficient  as  well  as  size  of  the  supersonic 
region  is  quantified. 

Tests  of  the  adequacy  of  the  foregoing  model 
to  simulate  realistic  transonic  USB  airfoil  flows 
have  been  inhibited  by  the  lack  of  suitable 
experimental  data.  Information  exists  only  'or 
highly  three-dimensional  configurations,  large 
thickness,  or  incidence  in  ranges  beyond  the 
validity  of  the  assumptions  of  small  disturbance 
theory.  Another  restriction  is  the  unavailability 
of  the  associated  geometric  data  and  'low  diagnos¬ 
tics  accompanying  the  tests.  The  results  of 
Yoshihara  and  his  coworkers  were  useful  in  this 
connection  and  allowed  us  to  compare  the  jet  flap 
specialized  version  of  the  US3  theory  in  Ref.  2. 
ror  the  simulations  described  in  this  paper,  tests 
performed  by  N.C.  Freeman  at  NPL  on  a  US3  modified 
5S  thick  RAE  102  airfoil  and  described  in  Ref.  3, 
appear  to  be  the  most  suitable  results  'or  com¬ 
parison  at  present.  Unfortunately,  the  angle  of 
attack  associated  with  the  NPL  data  is  5’  which 
is  marginal  for  the  application  of  a  small 
disturbance  model . 

Figure  13  indicates  comparisons  of  chordwise 
pressures  for  various  values  of  Cj.  Also  shown 
are  schlierens  indicating  the  associated  flow 
field  structure.  Turning  to  the  Cj*0  -esults 
(Part  (a)),  massive  shock  induced  separation  is 
indicated  and  is  apparently  initiated  at  the  down¬ 
stream  limb  of  the  lambda  snock  on  the  ^oper  sur¬ 
face.  This  is  reflected  in  the  classical  erosion 
of  the  suction  plateau  and  is  'esoonsible  'or  the 
indicated  disagreement  between  the  inviscid  compu¬ 
tational  results  and  the  data.  ror  these  tests, 
nominal  tangential  blowing  with  a  slot  height  of 
3.37',  of  the  chord  was  used.  The  slot  location  is 
15T  downstream  of  the  nose.  The  Macn  number  M 
immediately  above  the  slip  line  at  the  slot 
(point  A  in  r1g.  1)  is  approximately  1.29  'or 
both  Cj’s  indicated.  ror  C--*0.D17,  the  slot  Mach 
number'Mo  has  been  estimated  as  1.79  ano  for 
C<  *  0.34$,  m  *2.36. 


Comoan  son  between  t.neory  and  experiment  in 
Part  (b)  of  Fig.  13  indicates  reduced  discrepancies 
on  tbe  upper  surface  associated  with  the  limited 
separation.  In  »art  (c),  the  agreement  is  corres¬ 
pondingly  further  improved. 

To  achieve  adequate  realism,  it  is  important  to 
discuss  factors  responsible  for  the  disagreements. 
One  'eature  not  captured  by  the  USB  simulation  is 
the  pressure  3 a i k e  at  the  slot  location.  Based  on 
the  slot  site,  the  streamwise  scale  for  this 
phenomenon  is  at  least  an  order  of  magnitude 
greater  than  the  characteristic  wavelength  of  a 
Mach  diamond  pattern  in  the  wall  jet.  These  fluc¬ 
tuations  may  not  be  resolvable  with  conventional 
pressure  tap  arrangements  fpr  the  thin  slot 
employed  in  the  tests.  If  a  rough  model  of  a 
coflawing  inviscid  supersonic  wall  jet  over  a  flat 
plate  is  used  to  describe  the  flow  near  the  slot, 
the  approach  to  a  final  steady  state  may  be 
damped  oscillatory  or  monotone  depending  on 
wnether  the  reflection  coefficient  R  which  is 
given  by 


where 


Turning  now  to  the  discrepancy  of  the  values 
shown  on  the  rear  surface  (downstream  of  0.5c)  in 
Fig.  13c,  we  note  that  in  spite  of  the  obvious 
elimination  of  separation,  a  thick  viscous  wall  jet 
is  present.  Downstream  diffusion  will  affect  the 
application  of  the  Soence  relation  on  the  blown 
portion  as  well  as  the  shock  jump.  In  view  of  the 
wall  jet  thickness  shown  on  the  schlierens,  this 
factor  appears  to  be  more  significant  than  shock 
obliqueness  at  its  foot.  A  near  term  refinement 
is  being  implemented  employing  a  coupled  inviscid- 
viscous  model  using  second  order  boundary  layer 
corrections  to  the  Soence  boundary  conditions 
accounting  for  axial  gradients  of  the  displacement 
and  momentum  thickness.  Once  these  refinements  ire 
incorporated,  systematic  optimization  between 
separation  suppression,  wave  drag  minimization,  and 
supercirculation  control  will  be  possible.  It  is 
envisioned  that  the  oesign  techniques  contained  in 
Refs.  13-15  will  augment  this  capability  by  provid¬ 
ing  methods  to  modulate  shock  formation  in  concert 
with  the  blowing  effects. 

A.O  Panel  us  ions 


Asymptotic  and  computational  models  have  been 
used  to  ootain  the  flow  over  upper  surface  blown 
(USB)  airfoils  in  incompressible  and  transonic 
flow.  The  treatment  involves  a  detailed  analysis 
of  the  flow  in  the  jet.  The  analytical  and  compu¬ 
tational  results  indicate  that 


is  respectively  positive  or  negative. 

The  relaxation  length  L  to  achieve  the  down¬ 
stream  pressure  in  units  of  the  axit  height  is  of 
the  order  of  in  A  which  can  be  approximately 
5  to  50  in  the  present  case  depending  on  the 
accuracy  of  the  estimate  for  Ma.  Note  in  this 
connection  that 


R<0  for  1  ,  and  M 

3e  9 

A >0  for  . 

"e 

For  the  submerged  case,  A-l,  (Mo»M),  and  the 
Frandtl  periodic  pattern  is  obtained,  with  no 
radiation  of  energy  to  the  external  flow. 

These  'acts  suggest  that  one  factor  that  may 
be  '“sponsible  for  the  observed  spike  is  the 
internal  decay  process  in  the  jet.  If  transonic 
effects  and  wall  curvature  are  accounted  for,  the 
presence  of  "ballooning"  and  throats  in  the  jet 
may  also  be  contributory.  We  have  discussed  such 
phenomena  in  connection  with  submerged  transonic 
wall  jets  in  Ref.  10  and  have  reported  analogous 
results  'or  the  coflowing  case  in  Ref.  11.  Selec¬ 
tion  rules  in  terms  of  Ma  and  M  for  the  existence 
of  throats  in  the  jet  near  field  are  given  in 
Ref.  12  wnich  are  based  on  an  integral  form  of  the 
<arman  3uderley  equation.  A  rough  sketch  of  the 
wave  system  that  could  explain  the  soikes  in 
Figs.  13b  and  13c  is  shown  -n  cig.  U.  Yet 
anothe'-  onencmenon  that  would  nave  a  similar  wave 
pattern  would  be  a  slight  upward  mot’on  of  the  jet 
due  to  viscous  mixing  or  a  misalignment  with  the 
surface  tangent  at  point  A. 


•  In  the  thin  jet  small  deflection  approxima¬ 
tion,  the  pressure  jumps  associated 

with  the  Spence  theory  prevail  even  if  the 
flow  is  rotational  and  compressible. 

•  The  asymptotic  developments  pravioeo  allow 
further  systematic  refinements. 

•  Effects  associated  with  initial  skewness 
and  vortlcity  inaccessible  to  other 
theories  can  be  assessed. 

•  The  dividing  streamline  leaves  tangent  to 
the  upoer  sur'ace  in  incompressible  flow. 

•  Computational  resuits  obtained  for  transonic 
'JS3  configurations  indicate  significant 
enhancements  of  lifting  oressures  associated 
with  blowing. 

•  Comparisons  with  experiment  indicate  the 
need  for  ’-efinements  incorporating  wave 
interaction  phenomena  near  the  jet  exit 
as  well  as  viscous  interaction  processes 
in  the  downstream  portion  of  the  wall  jet. 
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